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1 Introduction 

Phase transitions in QCD are currently very actively researched. Although 
most of the attention is focused on the phase transition at finite baryon 
density and temperature and the expected critical point in the phase dia- 
gram [U [2], there are still challenges even at vanishing density, such as a 
determination of the order of the phase transition or the flavor dependence 
of the chiral phase transtion temperature [3]. Since lattice QCD simulations 
have not given final answers to all these questions so far, complementary non- 
perturbative approaches are indispensable to gain a better understanding of 
the non-perturbative phenomena. 

For our non-perturbative studj0 of the chiral phase boundary of QCD at 
finite temperature, we use the functional renormalization group (RG) 
[3 [8] , applied to a formulation of QCD in terms of microscopic degrees of free- 
dom given by quarks and gluons. At low temperature and momentum scales, 
QCD can be described well by effective field theories in terms of ordinary 
hadronic states. But a hadronic picture is eventually bound to fail at higher 
temperature and momentum scales, owing to asymptotic freedom. The recent 
discussion of a strongly interacting high-temperature phase of QCD suggests 
that a simple description of QCD around and above the phase transition tem- 
perature does not exist [9]. In this case, a first-principles description in terms 
of quarks and gluons is most promising to bridge wide ranges in parameter 
space. 

Here we discuss particular problems which are accessible in a formulation 
in terms of quark and gluons [lOi lllj : first, we present our results for the 
running of the gauge coupling at finite temperature. Second, we discuss the 
interplay between gluodynamics and (induced) quark dynamics. Finally, we 
give results for the chiral phase boundary in the plane of temperature and 
number of quark flavors, obtained from such an investigation of the quark- 
gluon dynamics at finite temperature. 

The functional RG yields a flow equation for the effective average ac- 
tion r k [7], 



1 Talk given at the 2006 ECT* School "Renormalization Group and Effective Field 
Theory Approaches to Many-Body Systems" , Trento, Italy. 



2 Jens Braun 



d t r k = X - STrd t i? fe (/f 5 + i? fe r\ t = In |, (1) 

where Jj. interpolates between the bare action r k= A = 5 and the full quantum 
effective action r = -Tfc = o; -T^ denotes the second functional derivative with 
respect to the fluctuating field. The regulator function R k specifics the details 
of the Wilsonian momentum-shell integrations, such that the flow of 7^. is 
dominated by fluctuations with momenta p ~ k 2 . 

It is impossible to study the flow of the most general effective action, 
consisting of all operators that are compatible with the symmetries of the 
theory. Therefore we have to truncate the action to a subset of operators, 
which is not necessarily finite. Nevertheless, such an approximation of the 
full theory can also describe reliably non-perturbative physics, provided the 
relevant degrees of freedom in the form of RG relevant operators are kept in 
the ansatz for the effective action. This is obviously the most problematic part 
since it requires a lot of physical insight to make the correct physical choice. A 
first but highly nontrivial check of a solution to the flow equation is provided 
by a stability analysis of its RG flow, since insufficient truncations generically 
exhibit IR instabilities of Landau-pole type. 

The IR stability of RG flows can be improved by adjusting the regulator 
to the spectral flow of [12l [13] [14]. Doing this, we integrate over shells 

(2) 

of eigenvalues of rather than ordinary canonical momentum shells. In a 
perturbative language, the use of the spectrally adjusted regulator allows for 
a resummation of a larger class of diagrams. Such an improvement has been 
successfully applied to the study of QCD at zero and finite temperature [T2l 
1131 111] and is the underlying technical ingredient for the results presented 
below. 

In the subsequent sections, we use the exponential regulator [7] of the form 

f 2 } ("2") (2 

#fc( r fc ) = r k /[exptTfc /Z k k 2 ) - 1], where Z k denotes the wave function 
renormalization of the fields. 



2 Running gauge coupling at finite temperature 

In this section, we present our study of the running strong coupling at finite 
temperature. While an investigation of the running of the coupling is inter- 
esting in its own right, we will also show in the subsequent sections that it 
represents a key ingredient of our study of the chiral phase boundary. 

Owing to strong coupling, we cannot expect that low-energy gluodynamics 
are reliably described by a small number of gluonic operators. On the contrary, 
infinitely many operators become RG relevant in the low-momentum regime 
and drive the running of the coupling. We truncate the space of possible 
action functionals to an infinite but still tractable set of operators. For gauge 
invariance, we apply the background-field formalism as developed in [151 116) 
and follow the strategy of [TTJ [12] for an approximate resolution of the gauge 
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Fig. 1. (a): Anomalous dimension 77 as a function of a/8% = g 2 /327v 2 for SU(2) 
Yang-Mills theory at vanishing temperature. The anomalous dimension is defined 
via the background-field wave-function renormalization by 77 = —\ndtW^\ It is 
related to the flow of the coupling as dtg\ = r\g\. The arrow indicates the RG 
flow from the initial condition at the UV scale A (left dot) to the IR fixed point 
(right dot). Owing to asymtotic freedom, there exists also a fixed point for a = 0. 
(b): Running SU(3) Yang-Mills coupling atYm{k,T) as a function of k for tempera- 
tures T = 0, 100, 300 MeV, compared to the one-loop running for vanishing temper- 
ature. 



constraints [18] . For an introduction to the background- field formalism and 
to background-field flow equations for gauge theories, we refer the reader to 
the lecture notes of H. Gies [19]. Apart from standard gauge-fixing and ghost 
terms, our truncation of the gauge sector consists of an infinite set of operators 
given by powers of the Yang-Mills Lagrangian, 

r k [A, i>, ift] = J W k {9) + If [A, 4>, H 9= \F a ^F« v , (2) 

whereas the quark contributions are summarized in and will be discussed 
in Sec. [3] For the moment, we restrict our study to the gluonic subspace. Ex- 
panding the function W k {0) = W]p6 + jWf V + ^W^ } 9 3 the expansion 
coefficients span an infinite set of generalized couplings. The background- field 
method [20] provides us with a non-perturbative definition of the running 
coupling gk in terms of the background-field wave function renormalization 

Our truncation represents a gradient expansion in the field strength, 
which includes arbitrarily high gluonic correlators projected onto their small- 
momentum limit and onto the color and Lorentz structure arising from powers 
of F 2 . A drawback of such gradient expansions is the appearance of an IR un- 
stable Nielsen-Olesen mode in the spectrum [3_T[. At finite temperature T, 
such a mode will be populated by thermal fluctuations, typically spoiling per- 
turbative computations, see e. g. Ref. [55]. Our RG approach allows us to 
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resolve this problem with the aid of a temperature-dependent regulator which 
removes the unphysical thermal population of the unstable mode. Thus, we 
obtain a strictly positive thermal fluctuation spectrum. For details of the im- 
plementation of such a regulator, we refer to 

Inserting the truncation ([2]) in the flow equation (JTJ) , we find that the run- 

(i) 

ning of the coupling is successively driven by all generalized couplings Wf. . 

(i) 

Keeping track of all contributions from the flows of the , we obtain a non- 
perturbative (3 g 2 function in terms of an infinite asymptotic but resummable 
series in powers of g 2 , 



^^9 t9 i=^a m (i,N c )^^. (3) 



The coefficients a m depend on the temperature T and the rank N c of the 
gauge group0. At vanishing temperature, the P g 2 function agrees well with 
perturbation theory for small coupling. For larger coupling, we find that the 
integral representation of Eq. ([3j has a second zero, corresponding to an IR 
attractive non-Gaufiian fixed point g 2 > 0, in agreement with the results of 



Ref. [12], see Fig. 1(a) For an explicit representation of the a m 's and details of 
the computation, we refer the reader to Ref. [TT] . Note that the appearance of 
an IR fixed point in Yang-Mills theories is a well-investigated phenomenon in 
the Landau gauge [25, 26, 27], in accordance with the Kugo-Ojima and Gribov- 
Zwanziger confinement scenarios [28] [29]. Moreover, it is also compatible with 
the existence of a mass gap [T3] in Yang-Mills theory. 

As initial condition for the coupling flow, we use the value of the coupling 
measured at the r mass scale [3D], a s = 0.322. For scales k ^> T, we find 
agreement with the perturbative running coupling at zero temperature, as one 
would naively expect. In the IR, the running is strongly modified: The coupling 
increases towards lower scales until it develops a maximum near k ~ T . Below, 



the coupling decreases according to a power law g 2 ~ k/T, see Fig. 1(b) The 
reason for this behavior can be understood within the RG framework: first, 
the hard gluonic modes decouple from the RG flow at the scale k ~ T. At this 
point, the wavelength of fluctuations with momenta p 2 < T 2 is larger than 
the extent of the compactified Euclidean time direction. Hence these modes 
become effectively 3-dimcnsional and their limiting behavior is governed by 
the spatial 3d Yang-Mills theory. However, the decoupling of the hard modes 
alone cannot explain the decrease of the coupling for scales k < T. The second 
ingredient which is needed is the existence of a non-Gaufiian IR fixed point 
also in the reduced 3-dimcnsional theory. A straightforward matching between 



Odd powers in the coupling g which arise in high-temperature small-coupling 
expansions [23 1 124] are not reproduced in our calculation. This is due to the 
fact that we have not yet included non-local operators in our truncation. In any 
case, we do not expect that these operators are quantitatively important for the 
determination of the chiral phase boundary since, as we will discuss below, chiral 
symmetry breaking sets in at scales k <T. 
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the 4c? and 3d coupling reveals that the observed power law for the Ad coupling 
is a direct consequence of the strong-coupling IR behavior in the 3d theory, 
g 2 (k <C T) w g 2 d r k/T. Note that this asymptotic behavior can be deduced 
analytically from the the integral representation of Eq. ([3]), sec Ref. 
Again, the observed IR behavior at finite temperature is in accordance with 
recent results in the Landau gauge [31]. 

3 Chiral Phase Boundary of QCD 

For a study of chiral symmetry breaking at finite temperature, we extend our 
truncation @ in two directions. To make this more transparent, we divide 
if in Eq. ^ into two terms, if = r^ _kin ' + rf ~ int \ With the first term, we 
include the quark contributions to all gluonic operators, as done in Ref. [32] 
for QED and in Ref. [TT] for QCD: 



where we confine ourselves to massless quarks. An inclusion of mass terms is 
straightforward and has been done in Refs. This kinetic term successively 
contributes to the RG flow of the coupling and accounts for the screening 
nature of fermionic fluctuations. 

The second term, r^~ lnt ', accounts for the quark self-interactions. In a 
consistent and systematic operator expansion, the lowest non-trivial choice 
for this term is given bjH 



We have classified the four-fcrmion interactions according to their color and 
flavor structure: 



(V-A) = (VnvV*) 2 + (^7m75^) 2 , (V+A) = (^ 7 ^) 2 " (VHV75^) 2 , 
(S-P) = (V^) 2 - (^ 7 5^) 2 , (V-A) ad j = (^» 2 + (V^75T» 2 . 



Color and flavor singlets are given in the first line. Fundamental color (i, j, . . . ) 
and flavor £, . . . ) indices are contracted pairwise. The operators with non- 
singlet color or flavor structure are given in the second line, where (^ x t(j^) 2 = 
■tp x ij)^^ij} x , etc., and the (T a )ij denote the generators of the gauge group in 
the fundamental representation. 

3 Qualitatively, we expect our truncation of the fermionic self-interactions to be reli- 
able for our purposes, since the feed-back of higher-order operators, e. g. ~ (ipip) 4 , 
is in general suppressed by the one-loop structure of the underlying flow equation. 





+A VA [2(V-A) ad j + (1/JV C )(V-A)] . 



(5) 
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Fig. 2. (a): Representation of the one-particle irreducible (1PI) graphs which are 
contained in the RG flow equations for the four-fermion interactions. The double- 
lines represent the fully dressed propagator and the solid black dots denote the 
insertion of dtRk in the loop. The four-fermion couplings Ai are defined in Eq. |[5j. 
Note that diagrams with the same topology but with the regulator insertion attached 
to the other internal lines are present in the RG flow as well, (b) : Sketch of a typical 
j3 function for the fermionic self-interactions A; at vanishing temperature for zero 
gauge coupling (upper solid curve), for finite gauge coupling g > (middle/blue solid 
curve) and for gauge couplings larger than the critical coupling g > g CT (lower/green 
solid curve). The effect of finite ratios T/k is illustrated for vanishing gauge coupling 
by the dashed (red) line. 



In the following, we study the four-fermion couplings in the point-like 
limit, Ai = A.;(|pi| <C k). This is a severe approximation in the chirally broken 
regime where mesons manifest themselves as momentum singularities of the 
four-fermion couplings. Nevertheless, our point-like truncation can be a rea- 
sonable approximation in the chirally symmetric regim^f). We stress that our 
truncation ([5|) for the fermionic self-interactions forms a complete basis, i. e. 
any other pointlike four-fermion interaction which is invariant under SU(iV c ) 
gauge symmetry and SU(iVf)L X SU(iVf)R flavor symmetry is reducible by 
means of Ficrz transformations. Terms accounting for instantons, as Ua(1)- 
violating interactions have been neglected as well, since they may become 
relevant only inside the %SB regime or for small Nf. 

Introducing the dimensionless couplings Ai = fc 2 Ai, the (3 functions for the 
four-fermion couplings can be written in the form 

d t Xi = 2A, - XjAjkXk - bjXjg 2 - c t g A . (6) 

Here we use a straightforward finite-temperature generalization of the flow 
equations for the Ai, which have been derived and analyzed in |34l 133] . The 

4 This has recently been quantitatively confirmed for the zero-temperature chiral 
phase transition in many-flavor QCD [33 j . 
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quantities A, b and c are dependent on the ratio T/k, where A is a matrix, 
and b and c are vectors in the space of A couplings; see Ref. [TTj for an 
explicit representation of the flow equations. The various terms arising on the 
RHS of the flow equation of the A^s can be related straightforwardly to one- 
particle irreducible Feynman- graphs, which are depicted in Fig. 2(a) As initial 
conditions for the four-fermion couplings, we use Ai — > for k — > A — > oo. 
This choice ensures that the Aj are generated solely by quark-gluon dynamics 
from first principles. This point is very important, as it is contrary to, e.g., 
the Nambu-Jona-Lasinio model, where the four-fermion couplings serve as 
independent input parameters. 

Our truncation provides a simple picture for the chiral dynamics illus- 
trated by Fig. |2(b)[ At vanishing gauge coupling, we find fixed points for the 
four-fermion couplings at A,; 7^ and \ = 0, where the latter ones are IR 
attractive. By increasing the gauge coupling g, the RG flow generates quark 
self-interactions of order A ~ g i 1 resulting in a shift of the fixed points. More- 
over, we observe that the four-fermion couplings approach fixed points A* in 
the IR, if the gauge coupling in the IR remains smaller than a critical value g cr . 
At these fixed points, QCD remains in the chirally invariant phase. Increasing 
the gauge coupling beyond the critical coupling g > g CT , the gauge-fluctuation 
induced A's become strong enough to contribute as relevant operators to the 
RG flow. This is reflected in a destabilization of the IR fixed points [34l 133] . 
In this case, the four-fermion couplings increase rapidly and approach a di- 
vergence at a finite scale k = fc x sB- Indeed, this strong increase indicates the 
formation of chiral quark condensates and therefore the onset of chiral sym- 
metry breaking. We recall the NJL model as an illustration: There, the mass 
parameter m 2 of the bosonic fields in the partially bosonized action is inversely 
proportional to the four-fermion coupling, A ~ 1/m 2 . In addition, we know 
from such NJL-type models that their effective action in the bosonized form 
corresponds to a Ginzburg-Landau effective potential for the order parameter. 
In its simplest form, the order parameter is given by the expectation value of 
a scalar field. Symmetry breaking is then reflected in a non-trivial minimum 
of the this potential. Consequently, the scale /c x sb at which the four-fermion 
couplings diverge is a good measure for the chiral symmetry breaking scale. 
At this scale, the effective potential for the order parameter becomes flat and 
starts to develop a nonzero vacuum expectation value. 

At this point, we have traced the question of the onset of chiral symme- 
try breaking back to the strength of the coupling g relative to the critical 
coupling g cr . At zero temperature, the critical coupling for three colors and 
not too many (masslcss) quark flavors is much smaller than the IR fixed point 
value of the coupling. Therefore QCD exhibits broken chiral symmetry at zero 
temperature. 

At finite temperature, the running of the gauge coupling is strongly mod- 
ified in the IR. Moreover, the critical coupling becomes larger for increas- 
ing ratios T/k. This reflects that the formation of quark condensates require 
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Fig. 3. The figures show the running QCD coupling a s {k,T) for Nf = 3 massless 
quark flavors and N c = 3 colors and the critical value of the running coupling 
a cr (k,T) as a function of k for two values of the temperature, T = 130 MeV (left 
panel) and T = 220 MeV (right panel). The existence of the (a s ,a cr ) intersection 
point at the scale fc cr in the left panel indicates that the quark dynamics can become 
critical. 



stronger interactions at finite temperatures, since the quarks become stiffcr 
due to the thermal masses. 

In Fig. [3j we show the running coupling a s and its critical value a cr for 
T = 130 MeV and T = 220 MeV as a function of the regulator scale k. The 
intersection point k CT between both gives the scale where the chiral quark dy- 
namics become critical. Above this scale, the system is in the chiral symmetric 
regime, below it quickly runs into the broken regime. The critical temperature 
can now be estimated using the lowest temperature for which no intersection 
point between a s and a CI exists. We find T cr w 186 MeV for Nf — 2 and 
T a ?s 161 MeV for Nf — 3 massless quark flavors in good agreement with lat- 
tice simulations [3] ■ We stress that no other parameter except for the running 
coupling at the r mass scale has been used as input for the calculation of the 
critical temperatures. 

Finally, we discuss the chiral phase boundary in the plane of the tempera- 
ture and number Nf of massless quark flavors computed within our approach, 



see Fig. 4(a) For small Nf, we observe an almost linear decrease of the critical 
temperature for increasing Nf with a slope of AT a = T(Nf) — T(Nf + 1) f=a 
25 MeV. Additionally we find the existence of a critical number of quark 
flavors, iV f cr = 12, above which no chiral phase transition occurs. Since 
iV f cr < Nf = ^-N c = 16.5, our study confirms the existence of a regime 
where QCD is chiral symmetric but still asymptotically free. Our value for N" 
is in agreement with other studies based on the 2-loop /3-function [351 GEE E2] • 
A remarkable feature of the phase boundary is its shape near iV f cr . It is a 
generic prediction of the IR fixed-point scenario and can be understood from 
analytic arguments [11] : For Nf ~ N", where fermionic screening of color 
charges keeps the coupling small, the coupling must almost reach its maximal 
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Fig. 4. (a): Chiral phase-transition temperature T cr versus the number of massless 
quark flavors Nf. The flattening at Nf > 10 is a consequence of the IR fixed-point 
structure, see text. The analytic prediction ([7} for the phase boundary is depicted 
by the dashed (blue) line, (b): The figure shows the j3 function for Nf = 6,8,10 
massless quark flavors (from bottom to top). For illustration, we also show the linear 
approximation of the /3 function for iVf = 6 whose slope gives the corresponding 
critical exponent 0. 



value in order to drive the quark sector to criticality. This maximal value can 
be roughly estimated from the fixed point value g 2 of the coupling at zero 
temperature. In the fixed point regime, the (3 g 2 function is essentially given 
by the linear term of an expansion around the fixed point, /3 g 2 = —0 (g 2 — <?*), 
where denotes the universal critical exponent. For illustration, we show the 
(3 g 2 functions for Nf = 6, 8, 10 at zero temperature in Fig. |4(b)| Using the 
condition g 2 (k CI ) = g 2 r , we can relate the solution of the (3 g 2 function near the 
fixed point to the scale fc cr , where the coupling g 2 exceeds its critical value g 2 r . 
The fact that k cr sets the scale for the critical temperature T cr finally yields 
the following relation for the critical temperature near 7V f cr [TT] : 

T cr ~fc |iV f -iV f CT r*, (7) 

where should be evaluated at iV f cr and the scale fco is implicitly defined by a 
suitable definition for the (3 g 2 function. Relation ([7J is remarkable since it re- 
lates two universal quantities to each other: the shape of the phase boundary 
and the IR critical exponent 0. The good agreement between the full numer- 
ical result for the phase boundary and the analytic result Eq. ([7]) around N" 



is shown in Fig. 4(a) This allows us to conclude that the symmetry status of 
the system for a large number of quark flavors is governed by the fixed-point 
regime where dimensionful scales such as ^4qcd lose their importance. 
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4 Summary and Outlook 

We have shown that the functional RG allows for a systematic and consistent 
expansion of QCD. Our truncation of the effective action based on an oper- 
ator expansion seems to be promising for a study of QCD at zero and finite 
temperature, at least in the chirally symmetric regime. We have determined 
the chiral phase boundary of QCD in the plane of temperature and flavor 
number. Our quantitative results for small N{ are in agreement with lattice 
simulations for N{ — 2,3. For large Nf we observe a flattening of the phase 
boundary near iV f cr , owing to the IR fixed-point structure of QCD. Future 
extensions should include mesonic operators which can be treated with RG 
rebosonization techniques [55J 132] ■ We refer the reader to the lecture notes 
of H. Gies [IH] for an introduction to such techniques. An extension of our 
present approach in this direction would not only provide access to the bro- 
ken phase and mesonic properties, but also permit a study of the order of the 
phase transition. 

The author would like to thank H. Gies for his collaboration on the work 
presented here, and is grateful to him and to B. Klein for proofreading the 
manuscript. The author acknowledges support from the GSI Darmstadt, from 
the DFG under contract Gi 328/1-3 (Emmy-Noether program) and from a 
fellowship provided by the organizers of the ECT* school. 
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